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Systems of three and four quantum particles in the boundary-condition 
model are considered. The Faddeev-Yakubovsky approach is applied to con- 
struct the Fredholm-type integral equations for these systems in framework 
of the Potential theory. The boundary-value problems are formulated for the 
Faddeev-Yakubovsky components of wave functions. 



I. INTRODUCTION 

The report is concerned with a treatment of three- and four-body quantum systems with 
pair interactions described by the boundary conditions of various types. Models of a such 
kind including the model of interactions with hard core attract attention due to simplicity 
of description of particle interaction at small distances (see Refs. |J, @ for review). In these 
models, the repulsive part of the pair interaction is described by the boundary conditions 
for wave function which are set on a surface du in the (center of mass frame) two-body 
configuration space R n . Here, n is dimension of particles (usually n = 3 and du is a sphere 
in R 3 ). 

However, being really simple in the two-body case, the boundary condition model in the 
case of three and more particles gives a birth to certain mathematical difficulties untypical 
for few-body problems with smooth potentials. Thing is that one comes here to necessity 
to construct a scattering theory in exterior of noncompact surface formed by aggregate 
of cylinders du a x R^ -2 )™ being supports in the iV-body configuration space R^ -1 )™ for 
interactions in pair subsystems a. Standard equations of few-body scattering theory (see 
Ref. 0) were derived for non-singular interactions and are not adjusted to a work with the 
arising boundary value problems. There are few attempts || in physical literature to 
make an immediate regularization of these equations in three-body case using special limit 
procedures where at the beginning, one takes regular potentials with finite repulsive cores 
and then the cores are pulled to infinity. In their results, the papers [p, 0, are close 
to this approach, too. These attempts were not quite successful because resulting equations 
were not of Fredholm type. To make the problem unique soluble it is necessary to take into 
account additional considerations Q, [S]- 
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Another method is proposed in works |J- |TT| of S.P.Merkuriev and author. Following 
the traditional approach to treatment of the boundary value problems for elliptic operators, 
we use the Potential theory (see Ref. [|L2]| ) reducing the problems to studies of integral 
equations on the surface where boundary conditions are set on. As mentioned above, this 
surface in the case of A-particle system with N>3 is unbounded. Therefore, in a difference 
to the boundary value problems for compact surfaces, the equations of the Potential theory 
are not of Fredholm type. To transform them into those of Fredholm type we apply in 
three-body case, the method by L.D.Faddeev | 13[j . Namely we extract and inverse explicitly 



singular diagonal part of integral operator in the Potential equations. As a result we obtain 
the Faddeev-type equations for densities of simple and/or double layer (kind of densities 
depends on a type of boundary conditions) |§- |T| which are did of Fredholm type and 
make possible to ground the scattering problem JTT]. Note, in the recent works |14f . 



a new 



method to prove the completeness of wave operators in iV-body problems with arbitrary N 
has been developed. This method is based on the concept of locally conjugate operator and 



is extended on the hard-core A-body Hamiltonians [15 



Alongside with the integral form we use also a differential form of the Potential equa- 



tions. Here, a very important thing are the generalized potentials (quasipotentials) @, flT 
(one-dimensional variants of quasipotentials for the boundary-condintion model were con- 
structed for the first time in Refs. JTjJ, |fPT| ) allowing to reformulate boundary conditions 
in terms of singular distributions. These potentials are presented by linear combinations of 
the delta-function concentrated on the surface and its derivative with respect to the surface 
normal. The generalized potential method simplifies a scheme of derivation of the Faddeev- 
Yakubovsky-type integral (in framework of the Potential theory) as well as differential (in 
terms of quasipotentials) equations for components of resolvent and then for the wave oper- 
ators. In the work we formulate such equations in the cases of three and four particles. 
However the derivation may be spread also on the case of a system with arbitrary number of 
particles. As in ||, integral equations allow to study asymptotical boundary conditions for 
the Yakubovsky components of wave functions at large values of space variables. The formu- 
lations obtained of the boundary value problems for the Faddeev- Yakubovsky differential 
equations may be applied to study of concrete few-body systems in the boundary-condition 
model. Some results of computations of three-body scattering and bound-state energies are 
presented in Refs. 0, @. 



In the present work, we review the results from Refs. |J- |TT[ , |2fl and [|TI| reformulating 
them for the case of arbitrary dimension of particles n>2. 



II. NOTATIONS 

In the report, we restrict ourselves only to the cases of N= 3 and N= 4 particles. For 
the sake of simplicity we suppose all the particles to be spinless. 

For description of three-body system we shall use the standard reduced relative coordi- 
nates || x a ,y a , a = 1, 2, 3. For example in the case of a = 1 these coordinates are expressed 
through the radius-vectors e R n and masses m^ of particles by the formulae 



2m 2 m 3 



.1/2 



m 2 + m 3 



(r 2 -r 3 ), i/i 



2m 1 (m 2 + m 3 ) 
mi + m 2 + m 3 



1/2 



ri - 



m 2 r 2 + m 3 r 3 * 
m 2 + m 3 , 
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Relative coordinates are combined in 2n-vectors X = (x a ,y a ). A choice of coordinate pair 
fixes cartesian coordinate system in R 2n . 

In the case of four-body system, alongside with the index a denoting again a pair sub- 
system, we use also the index a for partition of the system in two subsystems. If the pair 
a belongs to one of these subsystems, we write a C a. As above, we introduce the relative 
coordinates x a ,y a , x a G R n , y a G R 2n . Here, y a is a set of relative coordinates of the pair 
a considered as a whole and two particles in the rest. A detailed description of relative 
coordinates for four-body system can be found in Ref . [0] . 

In the boundary-condition model, the configuration space O of A-body system includes 
the points X G R^^ 1 )™ satisfying conditions x a G u a for all indices a where u a , u> a C R n , 
is the domain outside with respect to a piece- wise smooth closed compact surface du a C R n . 
This duj a is a surface where the boundary conditions are set on in the problem of two particles 
belonging to pair a. 

By r a we denote the [(N — l)n — l]-dimensional cylinders in R( Af_1 ) n engendered by 
surfaces du a} T a = du a x R (Ar_2)n = {XeR {N ~ 1)n : X = (x a ,y a ), x a G du a }. 
Hamiltonian of ^-body system is defined in 1/2(0) by the expression 

Hf(X) = f-Ax + ^faW) f(X) 
on the functions / G W|(0) satisfying the Dirichlet conditions (hard-core model) 

/!«! = « (!) 

or conditions of the third type 

= 0, x a e du a , a = 1,2,3, (2) 

on the boundary 90 of the domain O. Smooth functions parameters of the model 

and are defined for x a G du a . Potentials v a (x a ) describe pair interactions of particles at 
x a G uj a and are supposed to be smooth quickly decreasing functions. 







-\~ T a ( X c 



III. FADDEEV EQUATIONS FOR DENSITY OF SIMPLE LAYER 



At the beginning, we consider the first approach 0, JT0| to the boundary-condition model 
based immediately on the Potential theory jl2| . We demonstrate it for the case of N= 3 and 
conditions ([[]) supposing also that v a = 0, a — 1,2,3. In this case the kernel R(X,X',z), 
X,X' G O, of resolvent R(z) = (H — z)^ 1 satisfies the identity (following from the Green 
formula) : 

R(X,X',z) = R (X,X',z)~ I da s R (X,S,z)—R(X,X'z), (3) 

Jan on s 
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with Rq(X, S, z), Rq(z) = (—Ax — z)^ 1 , the Green function^ of the Laplacian —Ax in 
L 2 (R 2 ™) and ns the normal to dVt directed in Q. It follows from the representation (|]) that 
Green function R(z) is explicitly expressed in terms of its normal derivative fx(Z, X'z) = 

-pf!—R(S, X', z). The surface integral J dasRo(X, S, z)fj,(S) = U (X, z) in (H) is the potential 



fl2|| in all the points S where surface dQ is smooth 

dU 1 

hm 



of simple layer |L2| with density /i. This potential is known to have finite normal derivatives 

urface dfl is smooth, 

an 

Here, the up (low) sign — (+) corresponds to inside (outside) limit, i.e. X <eQ (XGR 2n \fi). 

Differentiating Equation (||]) with respect to the normal ng and taking into account 
relationship ([J), we obtain the following simple- layer potential equation for n(S): 

+ / ^-R (S,S f ,z)fi(S')da s/ = ^-R (S 1 X',z). (5) 

2 J on s oris 

an 

In this equation, the variables X' and z are fixed parameters. For compact integration 



surfaces, the potential integral equations are known to be of Fredholm type [[12] . Equation 
(|) however, is not of the Fredholm type because the surface dQ is unbounded. 

Let us construct the Faddeev-type |L3| equations for the density fi. 

First, we shall introduce some new notations. Let be the part of cylinder T a belonging 
to dQ, T e a = T a dn. It is clear that dtt = \J a T%. We shall denote by r a the part of T a lying 
inside of dQ, T 1 = T a \ T%. Restriction \x a = /i| rc will be called the Faddeev component 

a 

of the density fi. It is convenient to consider these components as functions defined on the 
total cylinders r a . Then, by /z^ (//Jj we shall denote the part of fj, a (S) defined on r® (TJJ. 
A concrete definition of internal parts /J> l a {S) of the densities [x a will be given in the following 
paragraph. 

In this notation, Equation (||) may be rewritten as the system of three coupled equations, 

\pili a = P*V a R - P*V a R £ (6) 

A /3=1 

with P° (Pa) the operator of multiplication on the characteristic function of the set (rjj 
and V a Ro(z) the integral operator with kernel \ a R (S, X' , z) = Ro(S, X , z), S G r a , 

X' G R 2n . On the other hand, we use the equations © to define the internal parts fj} a of 
the densities fi a replacing P^ with P a . Total components fi a satisfy the equations 

1 3 
-/j, a = V a Ro — V Q i? ^2 Vp- 

1 0=i 



1 Remember that for —Ax in L2(R l/ ), the Green function Ro(z) is given explicitly, Rq(X,X', z) 

( l~\ {\fz\X — X'\) 
-j I ) — (^~ 2 )/ 2 — — ^ with ii^ 1 ) the Hankel function of the first type. 
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Further, following by the Faddeev method, we transfer diagonal terms with (3 = a to the left 
sides and inverse operators ^1 + Y a Ro appearing there. As a result we obtain the following 
equations JjlJ, Q for densities fi a : 

fi a = V a R a + \~Vp a n l a - V a R a /•$■ ( 7 ) 

Here, R a {z) stands for the Green function of three-body system with only hard-core inter- 
action in pair a. Operator Y a p Q has the kernel 

V aPa (S, S', z) = J ^-JL-R a (S, S", z)R%\S», S>, z)da s with the Green function for 

the Neumann problem inside of T a . Both kernels *V a R a and V a p a are explicitly expressed 
in terms of two-body subsystem a. 

The equations obtained have been studied by methods of the Potential theory [11]. It 



was shown in particular that after some iterations, the integral operator corresponding to the 
right-hand part of ([F]) may be present as a sum of compact operator and another operator 
with norm smaller than one. The latter is engendered by a neighborhood of the ribs of 
the surface dQ formed by intersection of cylinders T a . Therefore, the Fredholm alternative 
may be applied to these equations and properties of the density \i may be investigated. 
These properties being known, we study all the necessary properties of the Green function 
R(z) Pi . The further procedure [|ll|] of constructing the wave operators and studying 
their properties (completeness, orthogonality, asymptotics etc.) is quite analogous to that 
designed for three-body problems with smooth potentials ||. 

Equations similar to (|7|), are obtained and studied also in the case of the conditions (0) 



and non-zero pair potentials v a [1 I 



IV. FORMALISM OF GENERALIZED POTENTIALS 

Another approach to construct the integral as well as differential equations for com- 
ponents of resolvent (and wave functions) uses the formalism of generalized potentials ||, 
[TT|| . Note that equation (0) may be considered as the Lippmann-Schwinger equation with 

(quasi)potential V acting as Vf = Son-^f c - Here, dgnfi stands for generalized function (dis- 



tribution) called simple layer ]nj and f° r the limit values on Oil (taking from ) of the 

normal derivative f. Later, we shall use also notations -X- f 1 for similar limit values taken 
_ on J on J 

from R( 7V_1 ) n \ fi 5 and f e , f 1 for respective limit values on d£l of the function / itself. The 
generalized function Sgnfj, acts on in accordance with the rule (</?, Sgnfi) = J dn das^{S)p(S). 
Analogous notations will be used also in the case when surface dQ is replaced with cylinders 

r Q . 

8 

Let us introduce the two-body generalized potentials V a , V a f = ^r a ^/ e , and consider 
instead of H the new "operator" H, 

H/(X) = -A/pO+5> a /(X), (8) 

a 

with X running all the space R^ -1 )™. Thereby we spread the domain of H on functions 
defined outside as well as inside of the surface dfl. According to (f|) one has to suppose 
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these functions and their derivatives to be continuous right up to surfaceR T a , a = 1,2,3. 
However the breaks f c — f 1 ^ z and -g-f e — jfcf 1 must be allowed when r a is crossed. 
Set of these functions in W 2 2 (R ( - 7V ~ 1 - )n \ U r a ) will be denoted by V. Action of the Laplacian 



-A on V has to be understood in a sense of distributions, 



-A/ = -A / + 5> o 



'df or 



dn dn 



E 



d_ 

dn 



SrJf-f 



(9) 



where A stands for the usual Laplacian. With (^j) one can easily see that substitution of 
V a into the Schrodinger equation leads to the two-sided boundary conditions on r a , 



dn 



0. 



/1r 



(10) 
(11) 



This means [ IT]] that if the spectral parameter z is out of (discrete) spectrum of respective 
boundary value problems for domains engendered inside of R^ -1 )"- \ VL due to intersection 
of cylinders T a then the spectral problem = z^> for H is equivalent to that for H. 
Therefore, the conditions ([!]) may be replaced outside of this spectrum with the generalized 
potentials V a . 

In the same way one can consider the model with third- type conditions (0). Generalized 



potentials in this case are following |T1| : V a f = —Sr a ( T af c ) + -g^(^r a f e ) (see also [[J). 

It is easily to include in this formalism usual potentials v a . The only necessary thing is 
the replacement of V a in fl8|) with V a + v a . 



V. DIFFERENTIAL EQUATIONS FOR COMPONENTS 

In the formalism of generalized potentials, the components and differential equations for 
them are constructed in the same way as in the case of usual potentials 0. 

For example, the Faddeev components of the bound-state A-body wave function \l/ which 
are defined as U a = —Ro(z)[v a + V a ]\I/, satisfy the Faddeev differential equations 



-A x +v a + V a -z)U a 



(V a + V a ) £ Up. 



At X $.{JT/3, these equations turn into habitual form B 

P 

(-A x +v a - z)U a = -v a U P- 



(12) 



Presence of generalized potentials generates two-sided boundary conditions for components 
U a on cylinders T a . In the hard-core model (0) these conditions may be written as follows: 



2 Excluding the points belonging to the intersections U (r^ p| F a ) of cylinders T a , a = 1,2, 3. 
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y dn 



0. 



with a running all the numbers of pairs (see Refs. 
conditions in the case of the model @) read as 




(13) 



ID[ , [0 for details). Analogous 



(14) 



In the four-body problem, we introduce after Faddeev components, those of Yakubovsky, 
U aa = —G a (v a + V a ) E Up, where G a (z) is the Green function corresponding to only 

interaction v a + V a . Remember that U a = E U aa and ^ = J2U a . 



aDa 



The components U aa satisfy at X 1J Tp the Yakubovsky differential equations 

P 

(-A X + Va - z)U aa + V a U la = -V a ^2 U/3a- 

■yCa /3Ca 

Boundary conditions for them may be written 01 in the form 



(15) 



dn 



\P B L =0 



(16) 



in the hard-core model ([!]) and 



dn 



(17) 



in the case of conditions (Q). Here, (a, a) runs all the chains of partitions and 

^aa = EU ia + E E ^6- 

/3Ca 

Components of the scattering wave functions in the boundary-condition model satisfy 
the same equations ( |I2"D and ( |T5D and boundary conditions ( ]T3] , ^6] ) or ( |T^| , |TT[ ) . For the scat- 
tering processes, asymptotical boundary conditions as X — > oo for the Faddeev components 
(N=3) and for the Yakubovsky components (N=4) are similar to those in the case of usual 
potentials || , |[8| . With these conditions, the boundary value problems ( [T^JTB] ) and (|15|,|T6|) 
(or ([12| ) [14]) and (p~5| , p~T|) ) become uniquely soluble (at energies lying out of the discrete spec- 
trum of the respective external problem for dfl and some internal problems for the domains 
formed by intersection of cylinders r a ). Results of concrete three-body computations of 
(2 — > 2,3) scattering processes and bound-state energies may be found in Refs. JTO, 
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